Ternary Social Networks: Dynamic Balance and Self-Organized Criticality 
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Antal et al. [Phys. Rev. E 72, 036121 (2005)] have studied the balance dynamics on the 
social networks. In this paper, based on the model proposed by Antal et al, we improve it and 
generalize the binary social networks to the ternary social networks. When the social networks get 
dynamically balanced, we obtain the distributions of each relation and the time needed for dynamic 
balance. Besides, we study the self-organized criticality on the ternary social networks based on our 
model. For the ternary social networks evolving to the sensitive state, any small disturbance may 
result in an avalanche. The occurrence of the avalanche satisfies the power-law form both spatially 
and temporally. Numerical results verify our theoretical expectations. 

PACS numbers: 05.65. +b, 87.23.Ge, 89.75.-k 
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I. INTRODUCTION 

Antal et al. [l[ have studied the balance dynamics on 
the social networks based on the notion of social balance 
[1, [^. In the networks, each node is connected to all the 
others, representing each person knows all the others in 
the society. Each edge in the networks has two values, 
-1-1 and -1. If the edge is -1-1, it means the two persons 
are friendly with each other. If the edge is -1, it means 
the two persons are hostile towards each other. At every 
step, they choose a triangle randomly from the network. 
If the product of the three edges is -1-1, the triangle is 
stable. Otherwise, if the product is -1, the triangle is un- 
stable. For the stable triangle, it satisfies (i) the friend 
of my friend being my friend; (ii) the enemy of my friend 
being my enemy; (iii) the friend of my enemy being my 
enemy; and (iv) the enemy of my enemy being my friend. 
The unstable triangles always try to be stable, but the 
final state of the network depends on the edge flipping 
probability p, which is set manually, li p > | , the net- 
work will reach the state of " paradise" , with all relations 
being friendly. Several studies around the balance dy- 
namics have been carried out, including the studies of 
the university class of triad dynamics 0] and the satisfi- 
ability problem of computer science 0, etc. 

In the first half of this paper, at first based on Ref. 
[H, we map the edge relations to the node relations with 
some relaxation, so there are two opposing opinions in the 
network. We suppose for each triad relation in the social 
networks, it will change from stable to unstable or the re- 
verse, caused by the change of some persons' opinions in 
it. While for the change, the person's opinion depending 
both on the other opinions in the triad relation and the 
opinions all around him. Next, considering in the social 
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networks, there being two opposing opinions is extreme, 
we introduce the neutral opinion. So the values of each 
node are generalized to -1-1, and -1, and the metastable 
triad relations appear. We find the densities of each triad 
relation associated with the Hamiltonian and the geomet- 
rical temperature which is determined by the structure 
of the network without manual setting. Based on our 
model, we obtain the distributions of each triad relation 
and the time needed for dynamic balance both for the 
binary and ternary cases. 



In the latter half of this paper, we study the self- 
organized criticality (SOC) [6|-@] on the ternary social 
networks. The previous studies of the SOC on complex 
networks are mainly based on the BTW model [1, [9l4T3| . 
In this paper, we establish our model with some differ- 
ences. As we know, in the society, if one relation changes 
from one state into anther, it may affect other relations 
associating with it. The effect may be big or small, de- 
pending on many internal factors. So in this paper, we 
simplify this phenomenon, with an eye to the triad bal- 
ance dynamics, and make two modifications compared to 
Sec. nil We find under specific conditions the SOC to be 
observable. Besides, for the small- world network [T5l.[l6|. 
we find the way of construction have an influence on the 
occurrence of the avalanche. We analyze it theoretically. 
Numerical results verify our theoretical expectations. 



This paper is organized as follows. In Sec. HIl we map 
the edge relations to the node relations and generalize the 
binary networks to the ternary ones. For both the binary 
and the ternary networks, we obtain the distributions 
of each triad relation and the time needed for dynamic 
balance. In Sec. IIIIl we study the SOC on the regular 
network and the small-world network respectively, and 
find out the small-world effect on the occurrence of the 
avalanche. We end this paper with conclusions in Sec. 
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nodes from the selected triangle at random, then let the 
selected node do spin flipping. The spin flipping proba- 
bility is defined as 



(2) 





FIG. 1: a\ — a4 stand for each type of the triad relations, and 
the corresponding densities in addition. a\ and 04 are stable 
relations, while and a-^ are unstable relations. 



II. DYNAMIC BALANCE 
A. The Node Model 

First we define the node relations in analogy with Ref. 
[l| with some relaxation. We propose each node in the 
network has two values, +1 and -1, standing for each 
person's two opposite standpoints in the society, or two 
opposing opinions. For a triad relation, if all nodes have 
the same value, i.e., all persons have the same opinion, 
the relation is stable. If there are different opinions in 
a triad relation, conflicts are easily to occur, so the re- 
lation is unstable. The stable and unstable triad rela- 
tions are given in Fig. [T] In our model, the stable and 
unstable triad relations have the probabilities to change 
into the other, which is more realistic for the social net- 
works. As we know in the society, stable relations may 
become unstable because of conflicts and the surround- 
ings, and unstable relations may become stable because 
of the same interest and the surroundings. Considering 
each node has two values which have the probabilities to 
change into the other, we study the social networks in 
association with the spin systems. In the following con- 
text, we establish our model in association with the Ising 
model and the generalized Glauber dynamics [13, [iB] • 

When associated with the spin systems, each node in 
the network standing for a spin. We propose the Hamil- 
tonian of each triangle to be 



where iJ^ = —aia'^Oj + OjGk + cr^ffc) stands for the fl- 
nal Hamiltonian after the spin flipping, with <j[ the final 
value of the selected spin, and ctj , the values of the 
other spins unchanged. It is noted that the spin flipping 
probability depends on the final value, independent of 
the initial one [l^. In addition, Z — X](t'=-i-i _i e"^^'^ 
is the normalizing factor and /3 is a geometrical temper- 
ature determined by the structure of the network. Since 
/3 is not a real temperature, and in order to embody the 
majority principle, we propose 

= 1/1 ^ a,l{N -\)-o\\^ 1/|M - a% (3) 

where TV is the size of the network, and M — 
X^i i'^ii^)) is the average value of all spins. So the 
spin flipping probability depends both on the triad rela- 
tion and the surroundings 18, 19]. 

The numbers of positive and negative nodes in the net- 
work are as follows, 



N+ = 



N- 



(3ai + 2a2 + as) N 

= y (3ai + 2a2 + 03) , 
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C% (02 + 2a3 + 804) 



(^2 



N 

— (02 -I- 2a3 -I- 3a4) (4) 



where N is the number of total nodes and stands for 
the triangle density, satisfying '^i ~ 1- The densi- 

ties of each triangle attached to a positive node are 



N 



N+ 



3ai + 2a2 



as 



(5) 



where i = 1, 2, 3, and Ci is the number of positive nodes 
in ai. In a similar way, the densities of each triangle 
attached to a negative node are 



N_ 

3 ^J'-'J 



N- 



02 + 203 + 3a4 ' 



(6) 



where j — 2,3,4, with Dj the number of negtive nodes 



(1) 



where Ui , aj and au stand for the three spins of the tri- 
angle. In the binary case, <Ji,aj and cr^ can take either 
of the two values, -1-1 and -1. Besides, a > is a real 
parameter, standing for the strength of coupling. So the 
Hamiltonian is the form of the Ising model. The evolving 
rules are as follows. At every step we choose a triangle 
from the network randomly, and choose one of the three 



For each node, the spin flipping probabilities from one 
value to the other are given by Eq. ([2]), i.e.. 



, \m- 



g~ (l + M) {-<^j+'^i<^k-(yk) ^ g (1-M) (g'j+O'jg'fc+g'fc) 

(7) 

So the total flipping probabilities for each node from one 
value to the other are as follows. 
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p(l^-l) = 



,-a/(l+J\/) 



2a, 



-a/(l+M) 



ai 



03 



„3a/(l + M) 



g-a/(l+M) ^ g3a/(l-M) 3 g-a/(l+J\/) ^. g-a/(l-M) 3 g3Q/(l + M) ^ g-a/(l-M) ' 
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p{-l ^ l).V^p(-l^l) 
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02 



p3a/(l-J\/) 



2a3 



-q/(1-M) 



-Q/(1-Af) 



3 g3a/(l-Af) ^ g-a/(l+A/) ' 3 g-Q/(l-M) _|_ g-Q/(l + Af) ' g-Q/(l- Af) _|_ g3a/(l+Af ) ' 



(9) 




FIG. 2: A plot of Tjv with A'' in the log-log coordinate. 
M(0) = —1, A is the size of the network, and Tjv is the 
time needed for dynamic balance. The numerical results are 
obtained by taking the average of 5 x 10^ simulations on the 
networks with size from 2 X 10^ to 2 X 10''. 



FIG. 3: as and a^ represent the stable relations, ae the un- 
stable relation, while ag — aio the metastable relations. 



where ai to 04 stand for each triangle density, and Ci 
and Di the number of positive and negative nodes in 
respectively. When the network gets dynamically bal- 
anced, the densities of each triad relation are as follows, 

13 3 1 

ai = g, 02 = g, 03 = -, 04 = -. (10) 

We note the network reaches an explicit state, indepen- 
dent of the parameter a. And in the dynamically bal- 
anced state n+ = n~ , where 71+ and n~ stand for the 
densities of positive and negative nodes respectively. 

The time needed for dynamic balance is 

rjv-iv^("\ (11) 

where N is the size of the network, and C{a) > 0, being 
a function of a. The numerical results are given in Fig. 
[5] More details are given in Appendix I A 11 



B. The Generalized Node Model 



There being only two opposing opinions in the social 
networks is extreme. We propose the neutral opinion 
should exist, which means the person holding this opinion 
is indifferent to other opinions. So we propose each node 
should have three values, -|-1, and -1, standing for the 
three opinions. If there are the same opinion in the triad 
relation except the neutral opinion, the relation is stable. 
If there are different opinions in the triad relation except 
the neutral opinion, the relation is unstable. Otherwise 
the relation is metastable. The simplest way is by judging 
the Hamiltonian (see Eq. If i?A < 0, it is stable. 

If Ha > 0, it is unstable. If Ha = 0, it is metastable. 
There are four stable relations, three unstable relations, 
and three metastable relations in the network, as shown 
by Fig. [Hand Fig. 13 

Compared with Eq. Q, the numbers of positive, neu- 
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tral and negative nodes in the network are as follows, 

N 



TV" 



3 

N 

y 

N 



(3ai + 2a2 + 03 + 205 + + ag) , 
(05 + flg + a? + + 2aQ + 3aio) , 



N = — (02 + 203 + 3a4 + ae + 207 + ag) , (12) 

o 

where ai-aio stand for the densities of each triangle. The 
densities of each triangle attached to a positive node are 
as follows, 



(13) 



3ai + 2a2 + a3 + 205 +06 + 08 



where « = 1, 2, 3, 5, 6, 8, and Ei is the number of positive 
nodes in a,; . In a similar way, the densities of each triangle 
attached to a neutral node are as follows, 



05 + flg + fly + 2as + 209 + 3aio ' 



(14) 



where j = 5, 6, 7, 8, 9, 10, and Fj is the number of neutral 
nodes in aj . The densities of each triangle attached to a 
negative node are as follows. 



Gkik 



02 + 2a3 + 3a4 + ag + 2a7 + ag 



(15) 



where k — 2,3,4,6,7,9, with Gk being the the number 
of negtive nodes in a^. 

We study the ternary network in association with the 
Potts spin systems. Since the values of each spin are 
more than two, the spin flipping mechanism changes into 
the spin transition mechanism. For the spin transition 
mechanism, the spin takes one of the three values, +1, 
and -1, depending both on the triad relation and the sur- 
roundings p^. Il9|. The Hamiltonian, the spin transition 
probabilities and the geometrical temperature are illus- 
trated by Eqs. ([T]), © and © respectively. The total 
probabilities from one value to the others for each node 
are as follows. 



p(-l^O) = 



p(-l^l) 



15(0^-1) 
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3 g |A/ + i+Af g 3 
ag 1 2a7 



g IMl 
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f e +e i-M ei"i + e^+'" +e 1- 
1 , ag 



3 1 + e i+M+gi-M 3 l + ei+J'^+e i-'-f 9 

3a a 

02 ei-J>^ 2a3 e ^-'•i 



e 



04- 



3 g|i\/| -|- g l + M ^gl-J\/ 3 g |A/| ^g l + A/ +g 1-M g|AJ"| ^ g l + M -)- g 1-A/ 

ag graw 2a7 e^T=T" ag 
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e i+M 
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2a8 e^TTTu 
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07 



g l + Af 



3 g 1 + A.f -|- g |A/| _|_ g 1-A/ 3 gl + A/ ^g|A/| _|_ g l^Af 

2a9 eiTAT 



3 g l + A/ + 1 + gl-Aif 3 gl + Af -|- 1 -f g l-Af 3 
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07 



d g l + A/ -|- g |A/| _|_ g 1-A/ 3 gl + A/ ^g|J/| -|- g 1-A/ 

2a9 e~ flio 



3 g l+A/ + 1 + ei-'J" 3 gl+A/ + 1 + e 1^**" 3 
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g l + A/ ^ g |A/| g 1-A 

205 e~rn" 
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p(l ^0) = ai- 
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2a2 



0-3 
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g l + A/ -|_g|A/| ^gl-A/ 
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3 e l+A/ 4- 1 -I- g 1-A/ 3 gl+A/ +1 -I- g 1-A/ 9 
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(16) 



But Eq. (fT6)) is correct only for M ^ 0. When A/ — > 0, the transition probabilities are defined as (the explana- 



FIG. 4: The numerical results are obtained by taking the 
average of 20 simulations on the network with lO** nodes. Ini- 
tially M(0) = -1/2, t/(0) = 1/2, and the network takes a 
long enough time to get dynamically balanced. 



FIG. 5: A plot of Tjv with TV in the log-log coordinate. 
A/(0) = -1/2, y(0) = 1/2. iV is the size of the network. 
Tjv is the time needed for dynamic balance. The numerical 
results are obtained by taking the average of 5 x 10^ simula- 
tions on the networks with size from 2 X 10^ to 2 X 10*. 



tions are given in Appendix I A 1|) 

p(-l^O) = 
p(O^l) = 
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(17) 



Since the transition probabilities are not continued at 
Af = 0, we cannot give the analytical results of n+ , rt^ , n'^ 
vifith a, so we turn to the numerical method. The results 
are shown by Fig. 21 We note when a — > 0, n"*" = n~ = 
nP = 1/3. And when a — +00, nP = 0, which means 
the neutral opinion disappears, and the ternary network 
comes back to the binary network. 

When the network gets dynamically balanced, the dis- 
tributions of each triangle are as follows, 



Ci 



1+71 m 



(18) 



where n"*" ~ n~ = x and vP = y. Besides, i 
1, 2, 3, 10, and l,m,n are the respective numbers of 
positive, neutral, negative nodes in a^, along with Cimn 
the corresponding combinatorial number. For example 
^"201 — 3. 

For the time needed for dynamic balance, we expect it 
to satisfy 



N 



N 



D{a) 



(19) 



where N is the size of the network and D{a) > 0, being 
a function of a. 



The numerical results of the time needed for dynamic 
balance are given in Fig. [5l More details are given in 
Appendix lA 21 



III. SELF-ORGANIZED CRITICALITY 



As is well known, SOC [6|-|8| is studied on manjf com- 
plex systems, including on the random graph [9], the 
scale-free network 10-12] and the small-world network 
[l3| . etc, which are mainly based on the BTW model [6|. 
Besides, there are the studies of the SOC on complex net- 
works based on other models, which we refer the readers 
to Ref. lj| for details. In this paper, we study the SOC 
on the ternary social networks based on our model and 
the evolving rules, along with two modifications com- 
pared to Sec. mi The first modification is the struc- 
ture of the network. Because for a completely connected 
network, any disturbance is globe, with no propagation, 
so the occurrence of the avalanche will not follow the 
power-law form either spatially or temporally. Besides, 
compared to the completely connected network, the reg- 
ular network or the small- world network [l^, [I^ is closer 
to the real structure of the society. So we take the latter 
forms for study. The second modification is the judge- 
ment of a triad relation. In order to embody the majority 
principle, we take the j udgem ent as the gauge invariance 
or the gauge variance 



-22], which is explained as fol- 
lows. Under this judgement, the triad relation evolves 
from one state to the other by changing node's value. 

If the value of a node being -I- 1, it means the person 
shows the friendly face in the triad relation, and always 
tires to keep the relation stable. If the value being -1, it 
means the person shows the hostile face, and always tries 
to destroy the stability of the triad relation. While the 
value being 0, it means the person does not care whether 
the triad relation is stable or not. So we propose if the 
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friendly faces are in the majority, the triad relation is in 
one state, or the gauge invariance state. While the hos- 
tile faces are in the majority, the triad relation is in the 
other state, or the gauge variance state. The explicit ex- 
pressions of the gauge invariance and the gauge variance 
states [H-dl] are as follows. 

For each triad relation, we propose if 

CTl + 0^2 + 0-3 > or (Ti = (72 = ^3 = 0, (20) 

it satisfies the gauge invariance. If 

ai + (7-2 + (T3 < or (Ji ^ a2 ^ 0-3 = 0, (21) 

the triad relation does not satisfy the gauge invariance. 
Here cti , 172 and 173 are the values of the three nodes of the 
triangle. So ai, 02, 05, as, aio are gauge invariance states, 
while a3, a4, ag, a7, ag are gauge variance ones, as shown 
by Fig. [Hand Fig. 13 

We study the ternary network in association with the 
Potts spin systems. Since each node has three values, 
the way of spin changing is the spin transition [l^. Be- 
sides, we propose the spin transition depends both on 
itself and the surroundings [H, [l^. Considering the 
SOC is determined by the structure of the system, we 
set the geometric temperature as fi = ^ > 0, where 
T = I X]<i j> '^j/fT-i ~ fj'li with cr- the final value of i, and 
J2<ij> '^j/iT-i the average value of all neighboring spins 
of i. Based on these considerations, we give the spin 
transition probability of each node as follows, 

W{cr, ^ al) = le-/3^('^^s<.,,>^,)_ (22) 

In Eq. (EU, li{a[,Y.<^,,>^i) = -^.'E<.,>^j, where 
51/<i j> '^j the sum of values of all neighboring spins 
of «, and Z ~ X]ct'=-i +1 e^'^>^<*'J>'^^ is the normaliz- 
ing factor. In order to avoid the denominator of /3 being 
zero, we propose if the final value of the spin results in 
T = 0, it happens with probability zero, while the other 
values happen with probability one. In this way, what- 
ever the initial condition is, the network will evolve to 
the sensitive state. 

The evolving rules are as follows. We set all triad rela- 
tions in the network to satisfy gauge invariance initially. 
The simplest way is to set all spins to be zero, which 
stands for all persons in the society being neutral with 
each other initially. After a long enough time of evolu- 
tion, the network will reach the sensitive state. Then any 
small disturbance may result in an avalanche. 

For any small disturbance: 

(1) Choose a triangle at random. If it is gauge variance, 
nothing happens. If it is gauge invariance, we choose 
one of the three nodes at random, and let the node do 
spin transition according to Eq. ^F2^ . If after the spin 
transition, the triangle is still gauge invariance, nothing 
happens. Otherwise, we store all the gauge invariance 
triangles attached to the selected triangle, and goto step 
(2). 




FIG. 6: There are 9 x 10* nodes in the two-dimensional lattice 
with equal length and width. It is one of the simplest two- 
dimensional regular networks satisfying the triad relations. 

(2) Choose all the stored triangles successively at ran- 
dom. Because of the interaction, first we judge whether 
the triangle is gauge invariance. If it is gauge variance, 
nothing happens. Otherwise, let one of the three nodes 
do spin transition. If the final triangle is gauge invari- 
ance, nothing happens, otherwise we find out all the 
gauge invariance triangles attached to the selected tri- 
angle. 

(3) Find out all the gauge invariance triangles by the 
combined action of the stored triangles in step (2). If the 
number of gauge invariance triangles to store is nonzero, 
goto step (2), otherwise goto step (1). 

For each avalanche, we define the size as the number of 
triangles stored in step (2) changed from the gauge invari- 
ance state to the other state, and plus the one changed 
in step (1). Besides, we define the time of avalanche as 
follows. If the number of stored triangles in step (2) 
changed from the gauge invariance state to the other 
state is nonzero, the time is increased by one, and plus 
the one changed in step (1). 

At first, we study the SOC on the regular network. The 
structure of the regular network is given in Fig. |51 The 
size distribution and the time distribution are given in 
Fig. [7] (a) and[7](b) respectively. In the following content, 
we keep the number of nodes in both the regular network 
and the small- world network constant as 9 x 10^. Besides, 
the numerical results are obtained by taking the average 
of 100 simulations. For each simulation, we execute 9 x 
10^ time of disturbance. 

In Fig. [7] (a) and Fig. [7| (b), the diagrams are curved 
when s and t are small. Because we set all spins to be 
zero initially. It takes a long time for the network to 
reach the sensitive state. While the diagrams curved at 
the ends is caused by the finite size effect. 

Next, we study the SOC on the small- world network 
p^ . [Tg} . The structure of the small- world network is 
given in Fig. [51 

When a triangle changes from gauge invariance to 
gauge variance, the value of the selected node is de- 
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FIG. 7: Plots of P{s) versus s and P{t) versus t on the regular 
network, (a) s is the size of avalanche. P{s) is the distribution 
function, satisfying the power-law form, with exponent -1.76. 
(b) t is the time of avalanche. P{t) is the distribution function, 
satisfying the power-law form, with exponent -2.38. 



FIG. 9: A plot of P{s) versus s on the small-world network, q 
is the edge adding probability, s is the size of the avalanche. 
P(s) is the distribution function, (a) q — 0.05. The exponent 
of the power-law distribution is -1.91. (b) q = 0.3. The 
exponent of the power-law distribution is -2.79 




(b) 



N=9*10', q=0.3 



FIG. 8: From a two-dimensional lattice with equal length 
and width, for each node in the network, it sends out an 
edge with probability q, and attaches to another node at ran- 
dom, then attaches to one neighbor of the selected node at 
random. Multi-connection and self-connection are avoided. 
Thus, by adding edges, more triangles are created. The man- 
ner of adding random edges is shown by nodes a,b,c. 

creased. Because of the small-world effect, it may result 
in the triangles far away compared to the regular net- 
work becoming of gauge variance. If the triangles far 
away compared to the regular network become of gauge 
variance, it will block the spread. So we expect the size of 
avalanche to decrease, and the exponent of the power-law 
distribution to decrease accordingly. Based on these anal- 
ysis, we set the edge adding probability q of the small- 
world from q = 0.05 to q = 0.3, and observe the exponent 
of the size distribution decrease from -1.91 to -2.79, as 
given in Fig. [9] (a) and Fig. |9] (b) respectively. 

As is well known, for the small- world network [isl. [l6j . 



FIG. 10: A plot of P{t) versus t on the small-world network. 
q is the edge adding probability, t is the time of avalanche. 
P{t) is the distribution function, (a) q — 0.05. The exponent 
of the power-law distribution is -2.65. (b) q = 0.3. 



when the edge adding probability q is big enough, the 
structure of the network becomes random. We note the 
diagrams are curved in Fig. [3] (a) and Fig. 1^1 (b) when s 
is small, but they take different forms. It means that the 
random networks are easier for the spread in the initial 
conditions. 

Then there is the time distribution of the avalanche 
on the small-world network, which is similar to the size 
distribution of the avalanche, as given in Fig. [10] (a) 
and Fig. [TO] (b) for comparison. But we note, when q 
is big enough, with the structure of the network becom- 
ing random, the time distribution no longer satisfies the 
power-law form. Besides, there is condensation at the 
end of the diagram, which we want to study in the later 
work. 

At last, for the construction of the small- world net- 
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10' 



10' 



N=9*10', q=0.05 
5=4.0, P(t)-t'"' 



(b) 



10' 



FIG. 11: A plot of -P(s) versus s on the small-world network 
considering the distance effect, q is the edge adding prob- 
ability, s is the size of avalanche. P{s) is the distribution 
function, (a) q = 0.05, S — 1.0, and the exponent of the 
power-law distribution is -1.85. (b) q = 0.05, 5 = 4.0, and the 
exponent of the power-law distribution is -1.69. 



FIG. 12: A plot of P{t) versus t on the small-world network 
considering the distance effect, q is the edge adding prob- 
ability, t is the time of avalanche. P{t) is the distribution 
function, (a) q = 0.05, S = 1.0, and the exponent of the 
power-law distribution is -2.56. (b) q = 0.05, S = 4.0, and the 
exponent of the power-law distribution is -2.22. 



work, when considering the factor of distance, the small- 
world effect on the triangles far away compared to the 
regular network is weakened, which may be easier for the 
spread. Besides, the number of local triangles are in- 
creased. So we expect the size of avalanche to increase, 
and the exponent of the size distribution to increase ac- 
cordingly. For the small- world network, when considering 
the distance effect, it is constructed as follows. From a 
regular network, as shown by Fig. |S1 for each node in 
the network, it sends out an edge with probability q and 
attaches to another node in the network. The probability 
of attachment between nodes i and j is 

P., = „ '%_s ■ (23) 

Here 6 > 0, and £ij stands for the smallest distance 
between i and j on the regular network. Then it at- 
taches to one neighbor of the selected node at random. 
Multi-connection and self-connection are forbidden. In 
this manner, more triangles are created. For the small- 
world network, we keep q = 0.05 constant, and set S — 1.0 
and S — 4.0 for comparison. For the size distribution of 
the avalanche, the exponent increases from -1.85 to -1.69, 
as shown by Fig. [TT](a) and Fig. [Tl](b). Meanwhile, for 
the time distribution of the avalanche, the exponent in- 
creases from -2.56 to -2.22, as shown by Fig. [T^] (a) and 

Fig.[n](b). 

In a word, the size distribution of the avalanche always 
satisfies the power-law form, no matter on the regular 
network, the small-world network or the random net- 
work. But for the time distribution of the avalanche, 
it deviates from the power-law form on the random net- 
work, which we want to do further study in the later 
work. 



IV. CONCLUSIONS 



We have studied the social networks based on our 
model and obtained meaningful results. First, in Sec. 
|TT]of this paper, based on Ref. fl'l, we map the edge re- 
lations to the node relations with some relaxation, along 
with introducing the neutral relation. So the social net- 
works change from the binary networks into the ternary 
ones. We suppose the change of each triad relation de- 
pending both on itself and the surroundings. So based on 
our model, when the social networks get dynamically bal- 
anced, we obtain the distributions of each triad relation 
and the time needed for dynamic balance. Besides, we 
find under the extreme condition a — > -l-oo, the ternary 
networks come back to the binary ones. Second, in Sec. 
mil of this paper, we study the SOC on the ternary social 
networks based on our model. Because we suppose, if one 
triad relation changes, it will affect other triad relations 
associating with it. The effect may be big or small, de- 
pending on many internal factors. While compared with 
Sec. |TT]of this paper, we make two modifications for our 
node model as follows. One is the structure of the social 
networks, which takes the form of regular or small-world 
network. The other is the judgement of a triad relation, 
which takes the form of the gauge invariance or the gauge 
variance. When the ternary social networks evolving to 
the sensitive state, any small disturbance may result in an 
avalanche. Since the number of triad relations changed in 
the avalanche being no scale preferred, it should satisfy 
the power-law distributions both spatially and tempo- 
rally. Third, for the small-world network, we find out 
the small-world effect on the occurrence of the avalanche 
both theoretically and numerically. 
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fll = ^,02 = -,03 = -,04 = -. (A4) 



Appendix A: Distributions of each triad relation and 
the time needed for dynamic balance 

1. The node model 

For the node model, each node has two spin values, +1 
and -1. We define one step of time as executing N spin 
flipping successively at random [1]. Because each spin 
flipping is linked up with {N — 1){N — 2)/2 variations of 
the triangles, the variations of the densities of each triad 
relation with time are as follows, 



da I 
~dt 
da2 
~dt 

da3 



~dt 



- p(-l- 


1) "a. - 




-!)<, 




-!)< + 






-p(l- 




-p{-l 


1) "a. 










-p(l- 




-p{-l 


1) "as 


= p{l^ 




p{-l^ 





(Al) 



When the network gets dynamically balanced, we have 
dai da2 da^ da^ 



dt 



dt 



dt 



dt 



= 0, 



and 



(A2) 



(A3) 



We note the distributions of each triad relation indepen- 
dent of the parameter a, and n"*" = in the dynamically 
balanced state. 

Defining at any time the density of positive nodes to 
be p, we obtain 



M = n+ - = 2p - 1, 
p = (l + M)/2, p{<x) = 



2' 



(A5) 



where j is the number of positive nodes in ai and C| is 
the combinatorial number. We have the new equations 
as follows, 



dn^ 
"dT 
dn~ 

~dr 

dM 
~dt 



p{-l^l)-p{l 



-1) 



p(l^-l)-p(-l^l), 

dri^ dn'^ 
~dt dT 

2[p(-l^l)-p(l^-l)]. 



(A6) 



Substitute Eqs. ®, ® and JM]) into Eq. (jMI, we 
obtain 



dM 



l + M\ fl-M 



1 - M 



„3a/(l-M) 



g3a/(l-M) _|_ g-a/(l+M) 



l + A/\ fl-M 



-a/(l-M) 



g-a/(l-Af) g-a/(l+M) 



,-a/{l-M) 



g-a/(l-M) _|_ g3a/(l+M) 



1 + M 



-a/il + M) 



l + M\ /l-M 



e 



-q/(1 + M) 



g-a/(l+M) _|_ g3Q/(l-Af) 

1 + M\ fl-M^ ^ 



g-a/(l + M) _j_ g-a/(l-M) 

r" 



„3q/(1 + M) 



53q/(1 + M) _^_ g-a/(l-M) 



(A7) 



For Eq. ()A7[) . when M{t) ~ 0, after some calculations. For the dynamic balance of the network, we propose 

we obtain the approximate result. 



M{t) ^ t-^/^^°'\ 
where C{a) > 0, being a function of a. 



(A8) 



\MiTN)\^0 



(A9) 
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So the time needed for dynamic balance is 

rAr-iV^("^ (AlO) 

In the simulations, both for the binary and the ternary 
cases, if the final value a[ results in the denominator 
of P (see Eq. ^) being zero, we propose it happens 
with probability zero, while the other values happen with 
probability one (see Eq. (|17p for example). In this man- 
ner, no matter what Af(0) (or M{0) and p(0)) is, the 



da I 
~dt 
da2 
~dt 



dag 
~dt 

daw 
dt 



When the network gets dynamically balanced, it satisfies 



where i = 1, 2, 3, 10, along with the normalization 

10 

= 1. (A13) 

i=l 

Although Eqs. (IJTTI) are very hard to solve, from ap- 
pendix lA H and the symmetric spin transition mechanism. 



network will reach the dynamically balanced state. 



2. The generalized model 

When each node has three spin values, the variations 
of the densities of each triad relation with time are as 
follows. 
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(All) 



we expect n"*" ~ n~ and M(oo) ~ in the dynamically 
balanced state. Defining n"*" ~ n~ = x and n° = y in the 
dynamically balanced state, the densities of each type of 
the triangles are 

a, ~ G™„a;'+"y'", (AM) 

where Z,m and n are the numbers of positive, neutral 
and negative nodes in a, respectively. C'lmn is the corre- 
sponding combinatorial number. For example C201 — 3. 
Because near M = 0, there are fluctuations. The spin 
transition probabilities are (Af 7^ 0) 



= p (-1 ^ 1) + p (0 ^ 1) nl, - [p (1 ^ 0) + p (1 ^ -1)] n+ , 

= p(l ^ -l)n+ +p{-l ^ l)n-,+piO ^ -l)n%+p{0 ^ 1) n^g - [p (1 ^ 0) + p (1 ^ -l)]n+ 
-[p(-1^0)+p(-l^l)]<2, 

= p(1^0)?i+ +p(-1^0)n-,+p(l^-l)n+ +p(0^-l)n^Jio-[p(0^-l)+p(0^1)]n° 

-[p(-1^0)+p(-l^l)]n-9, 
= p (1 ^ 0) n+ + p (-1 ^ 0) n-g - [p (0 ^ -1) + p (0 ^ 1)] n^.^. 



p(-l^O) 
p(-l^l) 
p(0^ -1) 
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(A15) 
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and (M ~> 0) 



Pi-l- 


-^0) 


= 0, p{1^0) 


-0, 






^ 2x'^y + 2xy^ 


, 

^ 2 ' 




-1) 


= 2x^y + 2xy'^ 


^ 2 ' 


p(l^ 


-1) 


= 2x^ + 2x^y -f 


xy^ 
2 ' 


Pi-l- 


-^1) 


= 2x^ + 2x'^y^ 


(A16) 



The densities of each triangle attached to a positive 
node are 



cix'^ + 2c2X^ + 2c^xy + C4?/^ 
(2x + yf 



(A17) 



where i = 1, 2, 3, 5, 6, 8, and X]j=i = '^j = 1- That 
is, only one of the coefficient Cj is nonzero. In a similar 
way, the densities of each triangle attached to a neutral 
node are 



cix + 2c2X + 2c^xy + C4y 
(2x + yf 



(A18) 



where j = 5, 6, 7, 8, 9, 10. And the densities of each tri- 
angle attached to a negative node are 



Cix^ + 2c2a;^ + 2c'ixy + c^y"^ 
(2a; + 2/)' 



(A19) 



where fc = 2,3,4,6,7,9. 

Because the transition probabilities (see Eq. (fTO)) and 
Eq. ([T7| ) are not continued at M = 0, we cannot obtain 
the analytical results for the distributions of n~ and 
rfi with a, so we turn to the numerical method, as shown 
by Fig. ^ We note when a -t-oo, the neutral opinion 
disappears, and the ternary network comes back to the 
binary network. That can be explained as follows. In the 
triad relation, the more people try to keep it stable, i.e., 
to show the same opinion with the others, the less likely 
the neutral opinion will exist. 
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